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ABSTRACT 

Stellar bars are the most common non-axisymmctric structures in galaxies and their 
impact on the evolution of disc galaxies at all cosmological times can be significant. 
Classical theory predicts that stellar discs are stabilized against bar formation if em- 
bedded in massive spheroidal dark matter halos. However, dark matter halos have 
been shown to alleviate the growth of bars through resonant gravitational interaction. 
Still, it remains unclear why some galaxies are barred and some are not. In this study, 
we demonstrate that co-rotating dark matter halos with spin parameters in the range 
of ^ Adm ^ 0.07 - which are a definite prediction of modern cosmological models - 
promote the formation of bars and boxy bulges and therefore can play an important 
role in the formation of pseudobulges in a kinematically hot dark matter dominated 
disc galaxies. We find continuous trends for models with higher halo spins: bars form 
more rapidly, the forming slow bars are stronger, and the final bars arc longer. After 
2 Gyrs of evolution the amplitude of the bar mode in a model with Adm = 0.05 is a 
factor of ^ times higher, A2/A0 = 0.55, than in the non-rotating halo model. After 
5 Gyrs the bar is ^ 2.5 times longer. The origin of this trend is that more rapidly 
spinning (co-rotating) halos provide a larger fraction of trailing dark matter particles 
that lag behind the disc bar and help growing the bar by taking away its angular 
momentum by resonant interactions. A counter-rotating halos suppresses the forma- 
tion of a bar in our models. We discuss potential consequences for forming galaxies 
at high-redshift and present day low mass galaxies which have converted only a small 
fraction of their baryons into stars. 

Key words: galaxies: halos - galaxies: spiral - galaxies: structure - galaxies: kine- 
matics and dynamics - galaxies: evolution 



1 INTRODUCTION 



that 



Optical and near-infrared observations indicate 
nearly 2/3 of disc galaxies in the local universe are 



strongly barred JEskridge et al]l2000l: lGrosb0l et all 12004 
iMenendez-Delmestre et al.ll2007l : lBarazza et al.ll2008l '). How- 
ever, the fraction of barred disc galaxies might have sig- 
nificantly increased from high redshift until the present 
day. Wheres som e studies indicate that it is nearly con- 
stant up to z ~ 1 jjogee et all |2004| : lElmegreen et al.l 12004 
iMarinov 



ther form as a result of swing-amplification of grav- 
itationally u nstabl e m = 2 modes in a stellar 
disc (Toomic ,19641: IColdreich fc Tremaind [l979l: iToomrd 
198ll:ICombes fc Sanderslll98ll : ISellwood fc Wilkinsonlll993l : 
Polyachenkd 120131') . triggere d by galaxy interactions 



and mergers ([Noguchil 119871 : ISarnes fc Hernguistl Il99ll : 
iMiwa fc Noguchil 1 19981 ') or through the cooperation of or- 



JogeJuOOTI ) other studies indicate a decrease in 



bar fraction t owards higher reds hift to values as low as 20% 
at z ~ 0.84 (jSheth et al.l 120081 ). A common message from 
these observational studies is that bars can form in their 
host stellar discs, embedded in diverse environments, and 
survive through 7 — 8 Gyrs of look-back time. 

The formation mechanisms of a bar have been 
mostly studied with N-body simulations. They can ei- 
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bital streams which eventu ally lead to formation of a bar 
(jEarn fc Lvnden-Bellj 19961). Interaction wit h dark matter 
halo substructures ( Romano-Diaz et al.ll2008l ) as well as the 
resonant gravitational interaction with the surrounding dark 
matter 'particles' ha s been shown to prom ote bar formation 
(lAthanassoula 200a : lDubinski et aLll2009h . Without interac- 
tions or in otherwise isolated environments, the formation 
and evolution of bars in disc galaxies strongly depends on 
the properties of the surrounding dark matter halo. It has 
been known for a long time that a dynamically hot halo com- 
ponent is required to stabilize an otherwise gravitationally 
unstable rota ting self-gravitating d isc ag ainst bar formation. 
According to lOstriker fc Peebles! l| 19731 ). a stellar disc will 
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go bar unstable if the ratio of total rotational kinetic energy 
to the potential energy exceeds 0.14 ± 0.03. If the disc is 
embedded in a - non-rotating dark matter halo, the poten- 
tial energy increases a nd the disc is stabilized against bar 
for mation (IHohll 19761). T his finding was again confirmed 
by Efstathiou et al.l (|l982l ) and quantified in terms of ob- 
servable parameters like disc scale length and circular veloc- 
ity. However, a study by (|Athanassoulall2002l ') indicated that 
strong bars can form in cold discs embedded in massive dark 
matter halos casti ng doubts on the ge neral applicability of 
the above criteria (|Athanassoulall2008l ) . 

Traditionally, the dark matter halo was treated as a 
rigid component which does not respond to the dynam- 
ics of the disc stars. However, 'live' dark matter halos en- 
hance the growth of bars in model galaxies by resonant 
interactions of disc stars with halo 'particles' taking away 
the disc's angular momentum ([Dcbattista fc ScUwood 19981; 
Athanassoula fc Misiriotis |2002| : iHollev-Bockelmann et al.l 
20051: IWeinberg fc Katj I2OO7I : [Ceverino fc KlvpinI l2007l: 
Saha et al.ll2012l ). The exchange of angular momentum be- 



tween the disc and halo has been a subject of active research 
over the last two decades and our understanding of the bar 
instability has improved significantly. 

Most studies focused on the formation of bars in kine- 
matically cold discs (which easily form bars) embedded in 
close to spherical non-rotating dark matter halos. The evo- 
lution of systems, especially dynamically hot discs which 
might even be dark matter dominated were left unex- 
plored, mainly because such systems were not expected to 
form bars at all. However, such systems exist. Low Sur- 
face Brightness (LSB) galaxies are generally isolated, with- 
out nearby co mpanions and are d ominated by dark mat- 
ter at all radii (|de Blok et al.ll2001f ). The evolution of such 
galaxies is expected to depend entirely on internal pro- 
cesses an d it remains unclear how their discs become ba r 
unstable (JMatthews fc Gallagheilll997l : iPohlen et al.ll2003h . 
as ~ 30% of LSB g al axies reported in the sample of 
Matthews fc Gallagheil (119971 ) are barred. The study of 
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( 20131 ) indicates, that hot stellar discs inside a 
massive non-rotating halo are stable against the formation 
of strong bars; they can form only weak bars over several 
Gyrs. In this paper, we test the hypothesis that halo rota- 
tion - in accordance with standard cosmological models - can 
promote the formation of strong bars even in the presence 
of massive dark matter halos. 

One of the definite predictions from large scale 
dark matter simulations with ACDM cosmological 
models is that dark matter halos acquire a certain 
amo unt of spin due to the cosm i c shear and merg- 
ers (ISteinmetz fc BartelmannI 19951 : iBuUock et al.l I2OO 



i 



IVitvitska et al.ll2002l : ISpringel et al.ll2005l : iBett et ahlboOT ) 
and they are hardly non-rotating. Their spin distribution 
peaks at a value Adm ~ 0.035 with a tail towards halos 
with even higher Adm-values. It has been sugg ested that 
rotating halos could enhance the growth of a bar (|WeinbereJ 
[l985|) as these halos have a higher fraction of dark matter 
particles on prograde orbits which can absorb disc angular 
momentum even more efficiently. We conduct a systematic 
investigation of the effect of halo spin on the evolution 
of a hot stellar disc embedded in a massive dark matter 
halo. The effect is tested with a galaxy model which is 
expected to be stable (or weakly unstable) against bar 
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Figure 1. Circular velocity of disc stars (red dashed), bulge stars 
(blue dotted) and dark matter (green dot-dashed) of our fiducial 
initial condition galaxy. The galaxy disc dominated by dark mat- 
ter (or a spheroidal bulge component) at all radii. 




Figure 2. Radial profile of the Toomrc Q parameter of the stellar 
disc. The dashed black line shows the initial condition at T=0. 
The other lines show the profiles after 2 Gyrs evolution for model 
A-l-00 (black solid), B-|-03 (green), C-l-05 (blue), D-|-07 (red), and 
the model with the counter-rotating halo E-05 (cyan). At any 
radius and at any time the stellar disc is stable in terms of Toomre 
Q> 1. 



formation according to most of the conventional criteria of 
bar instability. We find that the spin of a dark matter halo 
has substantial effect on the bar instability and thereby on 
the galaxy evolution and carry out a detailed investigation 
using - cosmologically motivated - spinning dark matter 
halos. 

The paper is organized as follows. Section [2] describes 
the construction of initial galaxy models with different halo 
spins and the basics of the simulation set up. Section [3] out- 
lines different conventional criteria of bar instability. Bar 
evolution and the morphology of the disc is presented in 
section |3] Section 14.21 and section 14.31 discuss the effect of 
halo spin on the Ostriker-Peebles criteria and bulge forma- 
tion respectively. Section[5]is devoted to angular momentum 
transfer and density wakes in the halo. Section |S] describes 
the resonances in the system. The discussion and conclusions 
are presented in section [7] 
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Figure 3. Time evolution (from left to right at T = 0, 2, 3, 4, 5 Gyrs) of the stellar surface density profiles for model with varying 
dark matter halo spins. From top to bottom the spin of the co-rotating (with the stellar disc) dark matter halos increases (models 
A+00, B-l-03, C-l-05, D-l-07). The bottom row shows the evolution of a model but with a counter-rotating halo (E-05). Models with more 
rapidly spinning co-rotating dark halos become bar unstable earlier and develop longer bars. Counter-rotating halos do not promote bar 
formation. 



2 GALAXY MODELS WITH SPINNING HALO 

We construct 5 equilibrium models of galaxy with initially 
spinning dark matter h alos (hereafter, ISDM) usi ng the 
self-consistent method of iKuiiken fc Dubinskil ([1995J). Each 
galaxy model consists of a live disc, a dark matter halo, and a 



classical bulge. Flattened models of ISDMs are constructed 
by reversing the velocities of halo particles with negative 
angular momenta perpendicular to the disc plane, which re- 
main s a valid solution o f the coUisionless Boltzmann equa- 
tion (|Lvnden-Belll Il96 J : iHernguistl [l99^ : ISaha fc Gerhard! 
12013 ). In this way, we have constructed one model with a 
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Figure 4. Time evolution of the peak of the m = 2 Fourier 
amplitude of the surface density. Red dashed line indicates D+07, 
blue dotted C+05, green dash=dot B+03, black solid line A+00, 
and cyan dash-dot-dot E-05. There is a clear trend for models 
with more rapidly spinning co-rotating halos to develop stronger 
bars at earlier evolutionary phases (e.g., D-l-07). 



non-rotating halo (A+00), three models with co-rotating ha- 
los with spin parameters Adm ~ 0.03, 0.05, and 0.07 (B-f 03, 
C-l-05, D+07) and one model with a counter-rotating halo 
(E-05) with a spin parameter of Adm ~ 0.05, similar to model 
C+05 but the halo is rotating in the opposite direction of 
the disc. For the definition of th e spin parameter (A dm) of 
the dark matter halos we follow (JBuUock et al.ll200l| l. 



Ad 



1 Y.i "^^l^i X «» 



(1) 



where rrii, fl and Vi are the mass, radius and velocity of 
a dark matter particle; V^i,. = GMvir/Rvir denotes the 
circular velocity at the virial radius (-Rvir). The values of 
Adm are shown in Table [T] and reasonably cover the range 
of spin distributions of dark matter halos from cosmologi- 
cal simulations. An analysis of the Millennium simulation 
ISpringel et al.ll2005l ) indicates that the halo sp in parame- 
ter on ly weakly depends on halo mass and shape (jBett et al.l 
120071 ). Therefore, we changed only the halo spin parameter 
keeping its mass and shape unchanged to single-out the ef- 
fect of halo angular momentum on the disc evolution. 

The disc has an exponentially declining surface density 
with a scale-length Ra, a constant scale- height h^ and mass 
Md- The outer radius of the disc is truncated at 7.0Rd with 
a truncation width of 0.5Rd within which the stellar den- 
sity smoothly drops to zero. The live dark matter halo is 
modelled with a distribution function given by the lowered 
Evans model (^Evans 19931) which has a constant density core 
and hence produce a nearly flat rotation curve. The initial 
flattening of the halo is q = 0.8; the mass and tidal radius 
of the halo are given by M;, = 10.32Md and Rh — 36Rd 
respectively. The preexisting classical bulge is represented 
by a King distribution function (jKing 196y). The mass and 
tidal radius of the bulge are Mt — O.lSA/j and Rt = 3.2Rd 
respectively. For rel evant details on model c onstruction, the 
reader is referred to lSaha et all (|2010l . 120121 '). 

In terms of astronomical units, Md ~ 1.7 x 10^" A/© 
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Figure 5. Same as in Fig.|4]but for the m = A Fourier amplitude 
of the surface density. 



Table 1. Initial spin of the dark matter halos, values of Ostriker- 
Pecbles parameter and properties of bar at 5 Gyr. The unit of 
pattern speed Qb is in kms^^kpc^^. 



Models 


^dm 


top 


Rbar/Rd 


Rcr/Rd 


Hb 


A+00 


0.000 


0.039 


1.32 


3.28 


18.3 


B+03 


0.032 


0.053 


1.61 


3.91 


15.0 


C+05 


0.053 


0.078 


1.78 


4.72 


12.2 


D+07 


0.074 


0.115 


2.60 


4.76 


12.0 


E-05 


-0.053 


0.078 


1.07 


3.32 


18.0 



and Rd — 3.0 kpc, so that the time unit is 24.0 Myr. The 
total mass of the dark matter halo and the bulge is Mh = 
1.73 X 1O"M0 and Mb = 3.0 x 10^ Mq, respectively. We 
use a total of 3.5 x 10® particles to simulate each model; 
assigning 1.2 x 10® to the disc, 1.8 x 10® to the halo, and 0.5 x 
10® to the bulge. The softening lengths for the disc, bulge 
and halo are 9.0, 5.9 and 23.6 pc, respectively, following 
iMcMillan fc DehnenI (^2007^. All simulation s were performed 
using the Gadget code ( Springel et al.ll200H ) with an opening 
angle oi6 = 0.7 for the tree code performance and evolved all 
the models for about 5.5 Gyr. The half-mass rotation period 
(dynamical time-scale) at 2.5Rd is 237 Myr (37.75 Myr). At 
the end of the simulation the total energy was conserved to 
about ~ 0.03%. 

In Fig. [l] we show the circular velocity curve for the 
galaxy model constructed as described above. The model 
is dominated by dark matter (and a spheroidal bulge com- 
ponent) right from the central region, similar to s ome giant 
LSB galaxies (|de Blok et al.ll200l| - |Lelli et aLlbOld ). The ini- 
tial stellar disc is kinematically hot and stable against local 
gravitational instability as c an be seen from the initial ra- 
dial profile of the Toomre Q (|Toomrelll964l ) parameter (see 
Fig.©, defined as Q{R) = ar{R)K{R)/3.36GT.{R), where ar 
is the radial velocity dispersion, k is the epicyclic frequency 
and E is the surface density of stars. 
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Figure 6. Bar size as a function of halo spin. Open circles indicate 
the sizes of the bars at T = 3 Gyr and filled circles at T = 
5 Gyr. Open circle of a particular color refers to the same model 
as denoted by the filled circle of the same color. Bars in models 
with larger halo spin are larger at all times. 



3 CRITERIA FOR BAR INSTABILITY 

Many simulations have shown that a strong ba r forms eas- 
ily in a cold, rota tion dominat ed stellar di sc (jHohJ Il97ll : 
ISellwood fc Wilkins on 1993; Sa ha et al.ll2012l . and references 
therein). Based on the study of th e stability of Maclau- 
rin discs (JBinnev fc Tremaind Il987|) and N-bo dy simula- 
tions of galactic discs, lOstriker fc Peebles! (|l973) suggested 
that a stellar disc would become bar unstable if the ratio 
top = Tiot/IW^I exceeds a critical value of 0.14±0.03. In our 
simulations, we compute the Ostriker-Peebles criterion for 
each model in a slightly different way. Assuming the virial 
theorem applies to the isol ated galaxy, it can be shown that 
JBinnev fc: TremainjFigSTf ) 



iOP — 1/(2 + rrand/Tmcan)- 



Trand ^ 1^ j {"" ^ v)'^ p{x)(fx 



is the kinetic energy in random motions and 
1 



Tmcan ^ 1^ j "" p(x)d X 



(2) 



(3) 



(4) 



denotes the kinetic energy in the mean (primarily rotational) 
motion. To compute top for each model galaxy, we need to 
add the contribution to the total kinetic energy budget from 
all the sub-components in the model e.g., bulge (B), disc (D) 
and halo (H), so that 

(5) 



-i rand " — -'randjB \ -^ rand,D ~r -trand,H) 

and 

-t mean — -tmcan.B \ -t mcan.D \ -tmcan.H- 



(6) 



Since the initial bulge and disc parameters are identical 
for all the models in our simulation, the variation in the 
initial top arise due to the variation only in the halo spin. 
It is interesting to note that the total kinetic energy in a 
galaxy is independent of the sense in which particles are 
moving. So the total kinetic energy in a halo with a fixed 
spin parameter is the same, independent of whether it is 
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Figure 7. Time evolution of the Ostriker-Peebles stability pa- 
rameter, ioP: (section [Sj for models A-l-00 to E-05. According to 
this criterion all models should be stable against bar formation. 
However, this criterion accounts for the fact that the models with 
halo co-rotation are more prone to instabilities (e.g. higher values 
for top)- 



co-rotating or counter-rotating with respect to the disc; top 
is, however, sensitive to the net amount of halo rotation. It 
is worth noting that D-l-07 has a higher value than other 
models and is likely to be unstable according to Ostriker- 
Peebles criteria. We will return to this issue in a later section. 

A similar global criterion, suggested bv Efstathiou et al.l 

(|1982D , states that if adn = ^peak/VCMd/Tid < 1.1, a self- 
gravitating stellar disc would become bar unstable. Here, 
Vpcak refers to the maximum of the rotational velocity. Ba- 
sically, the dimensionless parameter Qein measures the ratio 
of the total mass within a given radius to that in the stellar 
disc, such that a massive halo would provide stability against 
bar formation. For a bare exponential disc without any bulge 
and dark halo, the value of Ocin = 0.63 (jEfstathiou et al.l 
Il982h and would be subject to vigorous bar instability. For 
our initial conditions Qdn = 1.3 indicating clearly that all 
the galaxy models would be stable against bar formation. 
Note that the value of Qdn is insensitive to the halo rota- 
tion. 

Acc ording to | Toomrg (Il98ll) and a number of other 
studies (|Sellwood fc Wilkinsonlll993l : iDubinski et~aLl |2009| . 
and references therein), bars in N-body simulations are 
formed through the swing-amplification mechanism. For this 
mechanism to work, the disc needs to have no inner Lind- 
blad resonance (ILR), allowing the feedback loop to be in 
action. It has been shown that to achieve strong swing am- 
plification of an m = 2 bar mode , a stellar disc requires 
to have Q < 2 and 1 < X2 < 3 (jjulian fc Toomrd Il966l : 
iToomrdligsIl '). Here A2(7?) = Rhc^ [R) / {A-kGT,{R)) is the 
swing-amplification paramete r for an m = 2 mode in th e 
tight winding approximation (jGoldreich fc Tremaindll979l V 
For the models considered here, the value of X2 = 6.4 at 
2.2Rd which marks the location of the peak of the disc's 
circular velocity. According to all classical stability criteria, 
the stellar discs in our simulation should be linearly stable 
against strong bar formation. 
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4 BAR EVOLUTION - MORPHOLOGY 

In Fig. |3l we show time evolution of the face-on surface den- 
sity maps (from left to right) of the initially axisymmetric 
stellar disc embedded in a massive dark matter halo with 
different spin parameters (from top to bottom, see Table [T}. 
The model with a non-rotating halo (A-l-00) remains nearly 
axisymmetric for about 2 Gyrs; however, as the spin of the 
halo increases, the disc starts forming a strong bar even ear- 
lier. It is clearly visible that the model with the strongest 
halo rotation (D-l-07) develops the strongest bar. The model 
with the counter-rotating halo (E-05) develops a bar in a 
similar way to the non-rotating model (A+00). 

We have quantified this finding in Fig. |4] which shows 
the time evolution of the peak of the m = 2 Fourier compo- 
nent of the surface density normalized by the axisymmetric 
(m = 0) component, A2/A0 for all the models. There is a 
clear trend that a more rapidly spinning co-rotating dark 
halo promotes the early formation of a strong bar. Already 
after 2 Gyrs, the amplitude of the stellar bar in the halo 
with Xdm = 0.074 is higher by a factor of 7 compared to 
the non-rotating halo case. Although at 5 Gyr, the values of 
A2 for all the models are nearly the same, the bar sizes are 
different (see section [4.11 and Table [l|. From Fig. O a sys- 
tematic increase in the bar length is noticeable as the spin 
of the co-rotating halo increases. A large-scale bar is formed 
in the model D-l-07 after 2 Gyr and by the end of the simu- 
lation at 5 Gyr the bar extends beyond ~ 8 kpc (th e peak of 
the r otation curve, see Fig. [ij . Swing amplification (JToomrd 
Il98ll ) might have played an important role in the formation 
and growth of such a bar as our preliminary check showed 
no ILR (Inner Lindblad Resonances) which is crucial for es- 
tablishing the feedback loop in the stellar disc. Note that in 
most models such bars do not trigger any two-armed spiral 
feature, except transient spirals; primarily because of the 
higher values of Toomre Q and X2 parameter of the disc. 
As a bar grows, the combined effect of th e bar and transient 
spira l s heat the stars th roughout the disc (jjenkins fc BinnevI 
I1990I : ISaha et al.ll2010l ) and as a result of that the Toomre 
Q of the disc rises to even higher values (see Fig. [2}. This, 
in turn, reduces the possibility of growing spirals through 
swing amplification. 

Interestingly, the growth of a stellar bar in our counter- 
rotating halo model (E-05) is significantly suppressed com- 
pared to the co-rotating model with the same spin (C-f05). 
Substantial differences between C-l-05 and E-05 are clearly 
visible in the early phases of evolution (around 2 Gyr, Fig. 
[3} . The evolution of E-05 is more similar to the model with 
the non-rotating halo (A-l-00). 

In addition to a bar mode, the stellar discs in these 
models also develop higher order modes for which we restrict 
our analysis only to the m = 4 component. In Fig. O we 
show the time evolution of the m — 4 Fourier component of 
the disc surface density. It is apparent from the figure that 
m = 4 mode is significant in model D-l-07, the one with the 
highest halo spin in our sample. Not only a bar but also such 
a m = 4 component would take part in taking away angular 
momentum from the stellar disc and thereby accelerate the 
overall transformation into a bar. The stellar discs of all 5 
models do not grow any significant m — 1 mode (or higher- 
order odd modes). 



4.1 Bar size and halo spin 

After ~ 3 Gyr, all the models have grown a bar in their 
stellar discs. The size of a bar increases as time progresses 
by loosing angular momentum primarily to the dark matter 
halo. We compute the bar size in our model galaxies to un- 
derstand its variation with the dark matter halo spin. There 
exists a number of w ays to measure the size of a bar in a 
galaxy (|Erwinll2005l ). In N-body simulations, the bar size 
measurements can be obtained from the radial profile of the 
m = 2 Fourier component {A2{r)) and its phase angle vari- 
ation with radius. We take the full width at half maximum 
(FWHM) of A2{r) as the first measurement of a bar size. 
Theoretically, the phase angle of a bar does not vary with ra- 
dius but in reality, they do. We consider the radius at which 
the phase angle varies by ~ 5° as the second measurement of 
the bar size. We take the mean of these two measurements 
to specify the bar size in our simulations. In Fig. [6] we show 
the bar size for the different models at two different epochs 
(after 3 Gyrs and 5 Gyrs). It is very clear that at a given 
time the bar has - at all epochs - significantly larger sizes 
for halos with larger spin parameters. The trend seems to 
be even non-linear. Note that in the case of a fast corotating 
halo, the bar size has increased by ~ 50% in 2 Gyrs. 



4.2 Halo spin and the Ostriker-Peebles criterion 

lOstriker fc Peeblea (|l973l ) showed how an otherwise bar un- 
stable stellar disc becomes more stable against bar formation 
as the mass of a surrounding dark halo increases. We show 
the time evol ution of tpp for our models in Fig. [7] (simi- 
lar to Fig. 6 of lOstriker fc Peebles! (|l973t )'). The initial value 
of top for the model A-l-00 suggests that the disc should 
be stable against bar formation. Increasing the value of the 
halo spin (keeping its mass constant) makes the system more 
bar unstable and model D-l-07 has fop = 0.115 and therefore 
should be marginally bar unstable. Indeed, this model devel- 
ops the strongest early bar in our simulations which can be 
interpreted in the sense that the spin of the co-rotating halo 
increases the net rotational motion in the galaxy as a whole. 
However, other models (B+03 and C+05) also quickly de- 
velop bars and from our analysis we find that a hot stellar 
disc can become bar unstable within a reasonable time if 
top > 0.03 which also happens to be the most probable spin 
parameter for a ACDM halo. However, in order to properly 
understand bar formation in galaxies, one needs to take into 
account the combined effect of the mass and spin (angular 
momentum) of the dark matter halo. 

4.3 Bulge formation 

Numerous N-body simulations have demonstrated that 
kinematically cold and thin stellar discs easily form bars that 
can transform into a centr al bulge- like structure, someti mes 
termed pseudo-bulges (see lKormendv fc Kennicuttl2004l . for 
a review). Boxy/peanut-shaped bulges are a class of bulges 
primarily formed v ia the buckling (or fire-hose) instability 
of a strong bar dCombes et aT] Il990l: iRaha et all Il99ll: 



Hernguist fc Weinberg 1992: lAthanassoula fc MisiriotisI 



20021: loebattista et al.ll2004l : iKormendv fc Kennicuttll20o'l 
Saha et all I2OIOI. |2013|). A detailed orbital analysis by 
Pfenniger fc Friedlil ( 199ll ) shows that the 2 : 1 vertical 
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Figure 8. Edge-on galaxy morphology (stellar surface density) 
after 5 Gyrs of evolution (right) compared to the initial conditions 
(left). Prom top to bottom the spin of the dark matter halo in- 
creases; except for the bottom panel (counter-rotating halo model 
E-05). Spinning co-rotating dark matter halos promote the for- 
mation of peanut-shapcd/boxy bulges (i.e. D-l-07). 



resonance plays an important role in the formation of 
such a boxy/peanut bulge. The time-scale to build such a 
bulge-like component in an isolated disc galaxy depends 
on various parameters such as the dark matter content, 
the value of Toomre Q in the disc, and might as well be 
influe nced by a pre-existing 'classical' bulge (|Saha et al.l 
120121 ). For kinematically cold discs {Q ~ 1), building a 
pseudo-bulge can take only one Gyr or less. In hot {Q > 2) 
and dark matter dominated discs, the time-scale to form a 
pseudo-bulge can be very long; in som e cases it might even 
take a Hubble time (|Saha et al.ll2013h . 

In Fig. |8] we compare the final {T — 5 Gyrs) edge-on 
stellar morphologies (projected stellar surface density) of all 
galaxy models to the initial conditions {T = 0). The models 
with a non-rotating halo (A-l-00) and a counter-rotating halo 
(E-05) show bulge structures similar to the initial conditions, 
except for some amount of t hickening due to heating by the 
weak bar (|Saha et al.l 120101 ). Although the stellar discs in 
both models (A-l-00 and E-05) form moderately strong bars 
at the end, they did not go through a buckling instability 
ph ase and have not de veloped a boxy/peanut-shaped bulge; 
see ISaha et al.l (|2013l ) for a similar case in a dark matter 
dominated galaxy model. On the other hand, all models with 
a co-rotating halo develop bars that are strong enough to 
develop a buckling instability and form peanut-shaped/boxy 
bulges. The effect is strongest for the model with the highest 
dark matter spin (D+07). As the galaxy models host pre- 




Figure 9. Time evolution of specific angular momentum of disc 
stars (upper panel), bulge stars (middle panel) and halo particles 
(bottom panel) for all 5 galaxy models. 1/^(0) denotes the disc 
specific angular momentum at T = 0. The red line denotes the 
model D-l-07, blue C-l-05, green B-l-03, black A-l-00, and cyan E- 
05. The mass of the halo is 57 times higher than that of the bulge. 
As expected the model developing the stronger bar transfers most 
angular momentum (D+07). 



existing classical bulges the models with halo spin develop 
composite bulges with a different origin and morphology. 



5 TRANSFER OF ANGULAR MOMENTUM 
AND DYNAMICAL FRICTION 

Once a bar starts to form from disc stars, it grows via 
non-linear gravitational interactions involving the exchange 
of angular momentum with the outer parts of the disc, 
the surrounding dark matter halo and a pre-existing 
classical bulge (if present). The angular momentum is 
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Figure 11. Time evolution (from left to right) of dark matter density wakes in the equatorial plane of the discs for models with varying 
halo rotation (from top to bottom). Red (blue) indicates over-dens (under-dens) regions with respect to the mean density at a given 
radius. The dashed line indicates the position angle of the stellar bar and the solid line follows the peak overdensity in the halo. The 
model with the highest halo spin (D-l-07) develops the strongest trailing density wakes in the halo (e.g., fourth row, second and third 
panel from the left). Trailing wakes result in angular momentum loss by dynamical friction. 
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Figure 10. Time evolution of the pattern speeds of the bars, Qb > 
for all models. In models with more rapidly spinning halos the bar 
pattern speed is lower at late times. 



tran sferred by the Kravitational torque produced by the 
bar (|Lynden-Bell fc Kalnaiall972l ). Thereafter, it has been 
emph as ized by se veral a uthors llTremaine fc Weinberel 



19841: IWeinberd Il985l: iHernauist fc Weinberg 
Debattista fc SellwoodI I2OO0I: [Weinberg fc Katz 



Athanassoulal |2002|. l2003l: ISellwood fc Debattistal booe 



1992 



2002 



Dubinski erZI l2009l : ISaha et all I2OI2I ) that orbital reso- 
nances between disc stars and the particles representing 
the dark matter halos play a significant role in the transfer 
of angular momentum between the bar and the rest of the 
galaxy, with dark matter halo particles taking away most 
of the angular momentum from disc stars. In section 15.11 
we discuss, in detail, the physical process that allows the 
dark halo to do so. Most previous studies have investigated 
the transfer of angular momentum to a non-rotating dark 
matter halo. Here we quantify the angular momentum 
exchange between a bar and a rotating dark matter halo in 
more detail. 

In Fig. [§1 we show the time evolution of the specific an- 
gular momentum of the disc, bulge, and dark matter halo 
for all models. As mentioned above, the exchange of an- 
gular momentum between the stars and dark matter halo 
occurs when a resonant condition In + nil — mils = be- 
tween the bar and halo orbit is satisfied, Jls being the bar 
pattern speed. For the bar mode (m = 2), the important 
resonances are in this context {I, n, m) — (—1, 2, 2) i.e., ILR; 
{l,n,m) = (0,2,2) i.e., corotation. Independent of halo ro- 
tation, all stellar discs loose angular momentum which is 
deposited in the stellar bulge and the dark matter halo. 
The rate of angular momentum transfer correlates with the 
spin of the dark matter halo (or the strength of the bar). 
The amount of angular momentum transferred is highest in 
mod el D+07, in acco rdance with the semi-analytic results 
from i Weinberg! (| 19851 ). According to his calculation, the rate 
of angular momentum loss from the bar monotonically in- 
creases as the fraction of prograde orbits in the halo increases 
from to 1. The reason for this is that prograde orbits con- 
tribute predominantly to the (— 1, 2, 2), (0, 2, 2) resonances 
which are, in general, stronger than ( — 1, — 2, 2), (0, — 2, 2) 
resonances for a given I. For the same reason, one would ex- 
pect less angular momentum transfer to a halo with a large 



fraction of retrograde orbits. In model E-05 with a counter- 
rotating halo, we see that the angular momentum loss from 
the stellar disc is nearly negligible up to ~ 3 Gyrs and only 
after 4 Gyrs it is comparable to the case of a non-rotating 
halo. Interestingly, the counter-rotating halo looses angular 
momentum between ~ 1.5 — 4 Gyrs of evolution. Only for 
model D-l-07 there is an indication that the angular momen- 
tum transfer to the halo might saturate. 

As a result of the angular momentum loss from disc 
stars, the pattern speed of the bars decreases (Fig. llOp . Mod- 
els with stronger bars (more rapidly spinning dark matter 
halos) show a more rapid slowdown of the bar. The bar in 
the model with the fastest rotating halo (D-f 07) slows down 
faster during the initial phase of evolution; the pattern speed 
of the bar in D-l-07 is nearly half of that in A-l-00 (non- 
rotating halo) at ~ 2 Gyr. Non-linear effects are also appar- 
ent from this figure. In model D+07, the pattern speed of 
the bar remains nearly unchanged from ~ 2.5 Gyr onwards 
although the halo keeps absorbing angular momentum (see 
Fig. [9]). Since the bar in this model keeps growing, we think 
most of the angular momentum loss is invested into changing 
the moment of inertia of the bar according to: 



dt dt dt 



(7) 



where Ib and Qb are the moment-of-inertia and pattern 
speed of the bar. Notice also that the angular momentum 
gain by the pre-existing classical bulge (in model D-f07) 
nearly saturates beyond ~ 2.5 Gyr - a fact that correlates 
well with bar pattern speed slow down beyond ~ 2.5 Gyr 
in the same model. Overall, we find that a co-rotating halo 
with a reasonable spin parameter is capable of absorbing a 
significant amount of angular momentum from the disc (due 
to the formation of a stronger bar). As a result the bar can 
be slowed down by factor of 2 within about a Gyr, i.e. about 
27 dynamical time scale computed at the half-mass radius 
after its formation. 



5.1 Density wakes and the dark bar 

A rotating self-gravitating bar inside a live dark matter 
halo will create a density wake in the halo density dis- 
tribution that in turn exerts a torque on the bar and 
slows it dow n. This phenomenon is known as dynam- 
ical friction (| Chandrasekharl 119431 : iTremaine fc Weinberg! 
[l98J) and has been studied, in great detail, by sev- 
eral authors, especially in the context of bar-halo friction 
('Weinber g!! 19851: iDebattista fc Sellwoodll2000l: I Athanassoulal 
[200a : .Scllwood fc Debattistall2006l : I Weinberg fc Katzll2007l l 



Previous investigations of the bar-halo friction using N-body 
simulations, however, considered only non-rotating halos. In 
the following, we study the nature of density wakes in mod- 
els with rotating halos in more detail. 

We decompose the density distribution of the dark 
matter halos in different models into spherical harmonics 
ph{r,9,ifi) = pim{r)Yim{6 , (fi) . The density component P22 
corresponding to Z = 2, m = 2 is the strongest, although 
the contributions from / = 4, m = 2 and I = 4, tti = 4 
are non-negligible, especially in models D-f 07 and C-l-05. In 
Fig. 111! we show the density distribution p22 projected onto 
the equatorial plane {9 = 90°) for each model at different 
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Figure 12. Surface density of stars of each model galaxy shown at the time when the bar has just formed. Formation times are given in 
the upper right corner of each panel. The bar forms first (after 0.8 Gyrs) in the model with the highest halo spin (D+07, left) and last 
(after 2.8 Gyrs) in the model with the counter-rotating halo (F-05, right). 
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Figure 13. EUipticity and position angles (PA) as a function of 
the semi-major axis of the bars at the time of their formation 
(see Fig. ll2l l. All bars have similar ellipticity profiles but different 
position angles and pattern speeds (see Fig. IIOI I. One sma unit 
corresponds to 0.04iljj. 



times. In model A-f 00 (top row), we see a tiny bar-like den- 
sity wake produced at about 3 Gyr. By the same time the 
density distribution of the counter-rotating halo model (E- 
05) remains nearly featureless without any density wakes . 
In contrast, strong bar-hke density wakes are produced in 
all models with co-rotating halos; with increasing intensity 
as the spin parameter increases. The nature of the density 
wakes can be followed clearly in model D-l-07 (fourth row 
in Fig. Ilip . In the initial phase (less than ~ 2 Gyr), dark 
matter particles corresponding to the I = 2,m — 2 density 



wakes form a spiral-hke structure which lags behind the bar 
in the stellar disc. The phase angle of the spiral-like struc- 
ture with respect to the stellar bar varies with radius. As 
time progresses, the density wakes get stronger and forms 
a bar-like feature in the inner part with a nearly constant 
phase out to ~ 2 disc scale lengths. At larger radii the wakes 
still have a spiral-like structure and lag behind the stellar 
bar. The swarm of dark matter particles in the spiral-like 
wakes lagging behind the stellar bar are the primary source 
of dynamical friction. They exert torques on the bar and 
slow it down. At around 3 Gyr the spiral-like structure in 
model D-f 07 disappears, leaving behind a bar-like structure 
which is almost in co-rotation with the stellar bar. This fea- 
ture has been called a dark bar o r halo bar in the lite rature 
(|Hollev-Bockelmann et"aLl l2005l : lAthanassouTal |2007| ) . This 
dark bar is a resonance driven feature and would be hard 
to distinguish from the stellar bar as it has nearly the same 
pattern speed. The same process happens in all other models 
with co-rotating halos, however less rapidly for halos with 
lower spin parameters. 



6 RESONANCES AND SLOW BARS 

To understand the nature and dynamics of the bars forming 
in models with spinning halos, we identify the birth time of 
the bars by inspecting the maximum values of A2/A0, the 
radial variation of the projected ellipticity and the position 
angle of the semi- major axis of the bars. As explained above, 
dark matter halos with high spin parameter promote the 
early formation of bars, which is quantified here. In Fig. 1121 
we show the face-on stellar surface density maps for all mod- 
els at the corresponding birth time of the bars which we have 
defined in the following way: a bar is 'born' when the peak el- 
lipticity (tmax) along the bar semi-major axis becomes larger 
than 0.3, and the the corresponding PA is almost constant. 
In Fig. 1131 we show the ellipticity and the position angle 
as a function of the semi- major axis (sma). The aforemen- 
tioned range of the peak ellipticity roughly corresponds to 
{A2/Ao)max ~ 0.2 (see Fig. |3]). At formation, all bars have 
very similar shapes but they have different pattern speeds. 
We compute the size of the nascent bar using only the ra- 
dial variation of its PA in the same way as explained in sec- 
tion 14.11 This way, we get t he largest siz e of the bar, called 
as Lbar in the literature (e.g. lErwinll2005l ) , in all models and 
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Table 2. Properties of nascent bars. The unit of pattern speed 
Cfl is in kms~^lcpc~^. 



Models 


Lbar/Rd 


ne 


Rcr/Rd 


Rcr/Ltar 


A+00 


0.85 


26.9 


2.20 


2.58 


B+03 


0.90 


27.0 


2.18 


2.42 


C+05 


0.72 


29.0 


2.00 


2.77 


D+07 


1.10 


29.5 


1.95 


1.77 


E-05 


0.80 


24.5 


2.50 


3.12 



the corresponding values are given in Table [2] We use the 
estimate of Lbar to get a lower limit on the ratio of the co- 
rotation radius to the radial extent of the bar, Rcr/Rtar, 
which decides whether these initial bars are slow or fast in 
nature. A bar is termed 'fast' if 1.0 < Rcr/ Rbar < 1.4 and 
'slow ' for which the ratio is greater than 1.4 (jAguerri et al.l 
I2OO3I ). 

Next, we calculate the locations of the resonances in 
the stellar disc to see whether spinning halos preferentially 
create slow bars. Observations and simula tions (wihout halo 
rotation) bv lChemin fc Hernandea (|2009r ) already indicated 
that systems dominated by dark matter can only devlop slow 
bars. Fig. [TJ] shows the location of corotation resonances 
when the bar has just formed in the disc and towards the 
end of the simulation. We have checked that there was no 
ILR in the stellar disc initially. From Fig. IIUI we know that 
the nascent bar in model D-l-07 rotates faster than all others 
(open red circle in Fig. I14p with CR at ~ 6 kpc. The slowest 
rotating bar in model E-05 (couter-rotating halo, open cyan 
circle) pushes CR to ~ 7.5 kpc. Comparing, the bar sizes and 
CR radii we find that 1.77 < Rcr/Rtar < 3.1 which is larger 
than 1.4 for all models, i.e. above the limit for fast bars. 
This property does not change with time. After 5 Gyrs of 
evolution the bars are in the range 1.83 < Rcr/Rtar < 3.1. 
The values for bar sizes, CR radii and pattern speeds after 
5 Gyrs are quoted in Table [l] 

We conclude that all bars in the galaxy models pre- 
sented here form as slow bars and they remain slow bars 
as the galaxy evolves. The transformation of a slow bar 
to fast is unlikely. Increasing the spin of a halo (to more 
unlikely values) could possibly reduce the ratio Rcr/Rtar 
close to 1.4 leading to the possibility of harboring a fast 
bar in a dark matter dominated galaxy. It is worthwhile to 
point out that bars in early- type galaxies are generally fast, 
while they cou ld cithe r be slow or fast in la te-type galaxies 
lAguerri et al.„200a : .Rautiainen et al.llJOOsI '). 



7 SUMMARY AND DISCUSSION 

In this paper we have demonstrated that spinning dark mat- 
ter halos promote the formation of stellar bars, the most 
common and most intriguing non-axisymmetric feature in 
galaxies. Even dark matter dominated hot stellar discs, 
which are expected to be stable agains t bar formation ac- 
cordi n g to traditional criteria (see e.g., lOstriker fc Peebles! 
19731: JEfstathiou et all 1 19821 . however, see 




Athanassoulal 



20081 ') , can form strong bars depending on the spin of their 



R (kpc) 

Figure 14. Location of corotation resonances in the stellar disc 
which has just formed a bar and at a later time, 5 Gyr. The for- 
ward hatched region denotes the range of pattern speeds amongst 
the nascent bars (at the time of their formation, see Fig. I12I I and 
the back-hatched region for the evolved bars at 5 Gyr. The ranges 
of pattern speeds are given in Table [T] and Table [21 The solid 
black and dashed green curves refer to f2 for D-l-07 at those two 
instances and dash-dot curve for the model A-|-00 at 5 Gyr. Open 
circles refer to the sizes and pattern speeds of the nascent bars 
while filled circles are those of the evolved bars. The colors are 
identical to previous figures. 



co-rotating host halos. Stellar discs embedded in faster spin- 
ning halos - we have tested halos in the range < Adm < 0.07 
- develop their bars up to 2 Gyrs earlier with a higher 
{m = 2) Fourier amplitude (up to A2/A0 ~ 0.7) and the final 
bars are long er (up to a fact or of 2.5). This effect has been 
suggested by (| Weinberg! 19851 ) and is caused by more efficient 
angular momentum transfer from the disc to the spinning 
dark matter halos due to resonant interactions between disc 
stars and co-rotating halo particles. We have demonstrated 
that dynamical friction between a bar and a halo is stronger 
when the halo is co-rotating with the disc leading to sig- 
nificant transfer of angular momentum to the dark matter 
particles which form a bar-like structure in the inner region. 
As a result, a higher local density of dark matter particles 
can be expected wherever the stellar bar is located. This 
effect is weakest in counter-rotating halos (we have tested 
one model) which thus appear to suppress the formation of 
bars. The presence of counter-rotating halos could plausibly 
explain why some of the observe d massive and cold stellar 
discs did not form any bar at all (|Sheth et al.ll2012l '). 

Not all bars in our mo dels form boxy /peanut-shaped 
(pseudo-) bulges in the end (|Kormendv fc Kennicuttll2004l : 
I Athanassoulal I2OO5I ) . Although the m = 2 peak amplitudes 
of all the bars are nearly the same after 5 Gyr, only the 
discs embedded in co-rotating halos with spin parameters 
Adm > 0.03 went through a buckling instability phase. Discs 
embedded in faster spinning halos develop stronger bars and 
appear more box-shaped in the end. 

The impact of halo spin on bar formation in 
the idealized model discs discussed here is unambigu- 
ous and might have important consequences for the 
formation and evolution disc galaxies in ACDM cos- 
mologies. As discussed in section [ij the rotation of 
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dark matter halos is a definite predicition from iarge 
scale dark matter simulations with ACDM cosmologi- 
cal models dS tcinmctz & BartelmanrJ Il99 5: Bull ock et al 



Springel et al.l i2005: Rett et al 



2001 : IVitvitska ct al. 2002, 
20071). As most of the stellar discs are e xpected to 
co-ro t ate with their host halos (see e.g., I Sales et al.l 
l2012l : lAumer fc White! 120131 ). the formation of bars might 
be promoted in general. It has been shown by a 
number of studies that systems with already massive 
and dyna mically cold stellar discs become bar unsta- 
ble e a sily (ICombes fc Sanders! 198ll: Vale nzucla fc Klvpin! 
20031: !Dubinski et al.! !2009l: !Sellwood fc D cbattistal !2009r 



Saha et al.! !2012l : !Athanassoula et al.l !2013l ) and a rotating 



halo might only accelerate the instability in those systems. 
More interesting is the effect on systems which tra- 
ditionally were not expected to be bar unstable but are 
observed to have bars. In particular, present day low- 
mass galaxies whic h are, in general, observed to be dark 
matter dominated JBarazza et al.! !200^: !Sheth et al.l !2008l : 



!Chemin fc Hernandez! !2009! : !Nair fc Abraham! !2010l ) and 
in the ACDM framework - are expected to have converted 
only a very small fraction of the available baryons into star s 
JYang et al.l!2012! : !Moster et a"Ll!2013! : !Behroozi et~aLl!2012h . 
Another class of galaxies for which the above effect might 
be very relevant are low surface brightness galaxies (LSB). 
They are dark matter dominated, their stellar discs are rel- 
atively hot at a very low level of star formation despite hav- 
ing a similar amount of gas as in their high surface bright- 
ness (HSB) counterparts. LSB galaxies are generally found 
in isolated environments which excludes the possibility of 
str ong tidal interaction yet bars are found in thes e galax- 



ies l!Matthews fc Gallaghei!!l997!:!Pohlen et al 
ous N-body simulations by (|Maver fc WadslevI 



2003h . Previ- 



2004l l showed 



that LSB stellar discs embedded in CDM halos are gener- 
ally stable against bar formation. Our study with spinning 
halos sheds new light on the prospect of bar formation in 
LSB galaxies. 

Another interesting (but definitely more speculative) 
regime for this process to become important is at high 
redshift (z > 2) where star forming galaxie s have hot 
disc components iForster Schreiber et al.!!2009r ) and a sig- 
nific ant fraction of the baryons is still in the form of 
gas JTacconi et al.!!2012l ). Here rotating halos might pro- 
mote the early and rapid formation of bars with important 
con sequences for insta bility driven bulge evolution mod- 
els l!Genzel et al.l!2008! l and, eventually, the formation and 
early growth of supermassive black holes ([Bower et al.l!2006l : 
!Fanidakis et al.l!2012! : fHirschmann et al.!!2012! V 
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